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Abstract 

A fundamental aspect in performance engineering of wireless networks is optimizing the set of 

links that can be concurrently activated to meet given signal-to-interference-and-noise ratio (SINR) 

C/3 . thresholds. The solution of this combinatorial problem is the key element in scheduling and cross- 

O 

layer resource management. Previous works on link activation assume single-user decoding receivers, 

■ that treat interference in the same way as noise. In this paper, we assume multiuser decoding 
■<!::j- , receivers, which can cancel strongly interfering signals. As a result, in contrast to classical spatial 

, reuse, links being close to each other are more likely to be active simultaneously. Our goal here 

is to deliver a comprehensive theoretical and numerical study on optimal link activation under this 
novel setup, in order to provide insight into the gains from adopting interference cancellation. We 
therefore consider the optimal problem setting of successive interference cancellation (SIC), as well 
as the simpler, yet instructive, case of parallel interference cancellation (PIC). We prove that both 

■ problems are NP-hard and develop compact integer linear programming formulations that enable 

■ us to approach the global optimum solutions. We provide an extensive numerical performance 
evaluation, indicating that for low to medium SINR thresholds the improvement is quite substantial, 
especially with SIC, whereas for high SINR thresholds the improvement diminishes and both 
schemes perform equally well. 
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I. Introduction 



In wireless networking, determining the sets of links that can be active simultaneously 
is a cornerstone optimization task of combinatorial nature. For a link to be active, a given 
signal-to-interference-and-noise ratio (SINR) threshold must be met at the receiver, according 
to the physical connectivity model [[T9l . Within this domain, previous analyses assume 
that the communication system employs single-user decoding (SUD) receivers that treat 
interference as additive noise. For interference-limited scenarios, it is very unlikely that 
all links can be active at the same time. Hence, it is necessary to construct transmission 
schedules that orthogonalize link transmissions along some dimension of freedom, such as 
time. The schedule is composed by link subsets, each of which is a feasible solution to the link 
activation (LA) problem. Thus, for scheduling, repeatedly solving the LA problem becomes 
the dominant computational task. Intuitively, with SUD, a solution to the LA problem consists 
in links being spatially separated, as they generate little interference to each other. Thus, 
scheduling amounts to optimal spatial reuse of the time resource. For this reason, scheduling 
is also referred to as spatial time-division multiple access (STDMA) [|28l . Optimal LA has 
attracted a considerable amount of attention. Problem complexity and solution approximations 
have been addressed in [|2l, BH, [[TSl . [[36l . A recent algorithmic advance is presented in [12]. 
Research on scheduling, which uses LA as the building block, is extensive; see, e.g, [7], 
[fm . [|29l and references therein. In addition to scheduling, LA is an integral part of more 
complicated resource management problems jointly addressing scheduling and other resource 
control aspects, such as rate adaptation and power control, as well as routing, in ad hoc and 
mesh networks; see, e.g., [El, El, Il26]| . 

In the general problem setting of LA, each link is associated with a nonnegative weight, and 
the objective is to maximize the total weight of the active links. The weights may be used to 
reflect utility values of the links or queue sizes If32]| . A different view of weights comes from 
the column generation, method proposed in O, |I71, which has become the standard solution 
algorithm for scheduling as well as for joint scheduling, power control, and rate adaptation 
[fTOl . The algorithm decomposes the problem to a master problem and a subproblem, both 
of which are much more tractable than the original. Solving the subproblem constructs a 
feasible LA set. In the subproblem, the links are associated with prices coming from the 
linear programming dual, corresponding to the weights of our LA problem. A special case of 
the weights is a vector of ones; in this case, the objective becomes to maximize the cardinality 
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of the LA set. 

All aforementioned previous works on optimal LA have assumed SUD, for which interfer- 
ence is regarded as additive noise. In this work, we examine the problem of optimal LA under 
a novel setup; namely when the receivers have multiuser decoding (MUD) capability [33]. 
Note that, unlike noise, interference contains encoded information and hence is a structured 
signal. This is exploited by MUD receivers to perform interference cancellation (IC). That is, 
the receivers, before decoding the signal of interest, first decode the interfering signals they 
are able to and remove them from the received signal. For IC to take place, a receiver acts 
as though it is the intended receiver of the interfering signal. Therefore, an interfering signal 
can be cancelled, i.e., decoded at the rate it was actually transmitted, only if it is received 
with enough power in relation to the other transmissions, including the receiver's original 
signal of interest. In other words, the "interference-to-other-signals-and-noise" ratio (which 
is an intuitive but non-rigorous term in this context), must meet the SINR threshold of the 
interfering signal. With MUD, the effective SINR of the signal of interest is higher than the 
original SINR, with SUD, since the denominator now only contains the sum of the residual, 
i.e., undecoded, interference plus noise. Clearly, with MUD, concurrent activation of strongly 
interfering links becomes more likely, enabling activation patterns that are counter-intuitive in 
the conventional STDMA setting. The focus of our investigation is on the potential of IC in 
boosting the performance of LA. Because LA is a key element in many resource management 
problems, the investigation opens up new perspectives of these problems as well. 

The topic of implementing MUD receivers in real systems has recently gained interest, 
particularly in the low SINR domain using low-complexity algorithms; see, e.g., [fT6ll . Tech- 
nically, implementing IC is not a trivial task. A fundamental assumption in MUD is that 
the receivers have information (codebooks and modulation schemes) of the transmissions to 
be cancelled. Furthermore, the transmitters need to be synchronized in time and frequency. 
Finally, the receivers must estimate, with sufficient accuracy, the channels between themselves 
and all transmitters whose signals are trying to decode. For our work, we assume that MUD 
is carried out without any significant performance impairments, and examine it as an enabler 
of going beyond the conventionally known performance limits in wireless networking. Hence, 
the results we provide effectively constitute upper bounds on what can be achievable, for the 
considered setup, in practice. 

The significance of introducing MUD and more specifically IC to wireless networking is 
motivated by the fundamental, i.e., information-theoretic, studies of the so-called interference 
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channel, which accurately models the physical-layer interactions of the transmissions on 
coupled links. The capacity region of the interference channel is a long-standing open 
problem, even for the two-link case, dating back to [[T]|, lfT3l . [l22l . Up to now, it is only 
known in a few special cases; see, e.g., lIH, ll30l for some recent contributions. Two basic 
findings, regarding optimal treatment of interference in the two-link case, can be summarized 
as follows. When the interference is very weak, it can simply be treated as additive noise. 
When the interference is strong enough, it may be decoded and subtracted off from the 
received signal, leaving an interference-free signal containing only the signal of interest plus 
thermal noise. 

From a physical-layer perspective, the simple two-link setting above corresponds to a 
received signal consisting of X = S + I + N, where S is the signal of interest, with received 
power Ps and encoded with rate Rs, I is the interference signal with received power Pj 
encoded with rate Rj, and is the receiver noise with power i]. Assuming Gaussian signaling 
and capacity achieving codes, the interference / is "strong enough" to be decoded, treating 
the signal of interest S as additive noise, precisely if 



where 7/ = 2^' — 1 denotes the SINR threshold for decoding the interference signal /. If 
condition ([T]) holds, i.e., the "interference-to-other-signal-and-noise" ratio is at least 7/, / can 
be decoded perfectly and subtracted off from X. Then, decoding the signal of interest S 
is possible, provided that the interference-free part of X has sufficient signal-to-noise ratio 



where 7s = 2^^ — 1 denotes the SINR threshold for decoding signal S. By contrast, if the 
interference is not sufficiently strong for ([T]) to hold, then it must be treated as additive noise. 
In such a case, decoding of signal S is possible only when 



This way of reasoning can be extended to more than one interfering signals. Towards this 
end, we examine the effect of IC in scenarios with potentially many links in transmission. 
Our study has a clear focus on performance engineering in wireless networking with arbitrary 
topologies. In consequence, a thorough study of the gain of IC to LA is highly motivated, 
in view of the pervasiveness of the LA problem in resource management of many types of 
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wireless networks. In the multi-link setup that we consider, the optimal scheme is to allow 
every receiver perform IC successively, i.e., in multiple stages. In every stage, the receiver 
decodes one interfering signal, removes it from the received signal, and continues as long as 
there exists an undecoded interfering link whose received signal is strong enough in relation to 
the sum of the residual interference, the signal of interest, and noise. This scheme is referred 
to as successive IC (SIC). From an optimization standpoint, modeling SIC mathematically 
is very challenging, because the order in which cancellations take place is of significance. 
Clearly, enumerating the potential cancellation orders will not scale at all. Thus compact 
formulations that are able to deliver the optimal order are essential, especially under the 
physical connectivity model, which quantifies interference accurately. 

Alternatively, a simplified IC scheme is to consider only the cancellations that can be 
performed concurrently, in a single stage. In this scheme, when determining the possibility 
for the cancellation of an interfering link, all remaining transmissions, no matter whether or 
not they are also being examined for cancellation, are regarded as interference. We refer to 
this scheme as parallel IC (PIC). It is easily realized that some of the cancellations in SIC 
may not be possible in PIC; thus one can expect that the gain of the latter is less than that of 
the former. A further restriction is to allow at most one cancellation per receiver. This scheme, 
which we refer to as single-link IC (SLIC), poses additional limit on the performance gain. 
However, it is the simplest scheme for practical implementation and frequently captures most 
of the performance gain due to IC. In comparison to SIC, PIC and SLIC are much easier to 
formulate mathematically, as ordering is not relevant. 

In EH, we evaluated the potential of SLIC in the related problem of SINR balancing. That 
is, we considered as input the number of active links, let the transmit powers be variables and 
looked for the maximum SINR level that can be guaranteed to all links. In [3], we exploited 
link rate adaptation to maximize the benefits of IC to aggregate system throughput. In parallel, 
another set of authors has made a relevant contribution in the context of IC [|27ll . They 
considered a SIC-enabled system and introduced a greedy algorithm to construct schedules 
of bounded length in ad-hoc networks with MUD capabilities. There though, the interference 
is modeled using the protocol-based model of conflict graphs Il23ll . which simplifies the 
impact of interference, in comparison to the more accurate physical connectivity model that 
we are using. 

The overall aim of our work is to deliver a comprehensive theoretical and numerical study 
on optimal link activation under this novel setup in order to provide insight into the gains 
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from adopting interference cancellation. This is achieved through the following contributions: 

• First, we introduce and formalize the optimization problems of LA in wireless networks 
with PIC and SIC, focusing on the latter most challenging case. 

• Second, we prove that these optimization problems are NP-hard. 

• Third, we develop ILP formulations that enable us to approach the global optimum 
for problem sizes of practical interest and thus provide an effective benchmark for the 
potential of IC on LA. 

• Fourth, we present an extensive numerical performance evaluation that introduces insight 
into the maximum attainable gains of adopting IC. 

We show that for some of the test scenarios the improvement is substantial. The results 
indicate that strong interference can indeed be taken as a great advantage in designing new 
notions for scheduling and cross-layer resource allocation in wireless networking with MUD 
capabilities. 

The remainder of the paper is organized as follows. In Section [III we introduce the notation 
and formalize the novel optimization problems. In Section [nil we prove the theoretical 
complexity results. In Section |IVl we propose a compact ILP formulation for the LA problem 
with PIC having quadratic size to the number of links. For the most challenging problem of 
LA with SIC, we devote two sections. In Section |Vl we treat SIC under a common SINR 
threshold. By exploiting the problem structure, we show that the order of cancellations can be 
conveniently modeled and derive an ILP formulation of quadratic size. Then, in Section |VIl 
we consider individual SINR thresolds. For this case, we give an ILP formulation of cubic 
size. In Section IVII[ we present and discuss simulation results evaluating the performance of 
all proposed IC schemes. In Section IVIIIi we give conclusions and outline perspectives. 

II. Definition of Link Activation vv'ith Interference Cancellation 

Consider a wireless system of K pairs of transmitters and receivers, forming K directed 
links. The discussions in the forthcoming sections can be easily generalized to a network 
where the nodes can act as both transmitters or receivers. Let IC = {1, . . . , K} denote the 
set of links. The gain of the channel between the transmitter of link m and the receiver of 
link k is denoted by Gmk, for any two m,k E IC. The noise power is denoted by i] and, 
for simplicity, is assumed equal at all receivers. The SINR threshold of link k is denoted by 
7jt. Each link k is associated with a predefined positive activation weight Wk, reflecting its 
utility value or queue size or dual price. If a link is activated, its transmit power is given 
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and denoted by pk, for k E IC. A LA set is said to be feasible if the SINR thresholds of the 
links in the set are met under simultaneous transmission. All versions of the LA problems 
we consider have the same input that we formalize below. 

Input: A link set /C with the following parameters: transmit powers pk, SINR thresholds 7^, 
and link weights Wk, VA; G /C, and gain values Gmk, Vm, k E IC. 

Consider first the LA problem with the conventional assumption of SUD, where the 
interference is treated as additive noise. This is the baseline version of the LA problem 
in our comparisons; its output is formulated as follows. 

Problem LA-SUD: Optimal link activation with single-user decoding. 
Output: An activation set A C IC, maximizing Yl ^fe ^i^d satisfying the conditions: 

WkEA. (4) 



Yl PmGmk + V 
m€A\{k} 



This classical version of the LA problem can be represented by means of an ILP formu- 
lation; see, e.g., |I71, [fT2ll . A set of binary variables x^, V/c E IC, is used to indicate whether 
or not each of the links is active. The activation set is hence A = {k E IC : Xk = I}. In 
order to ease comparisons to the formulations that are introduced later, we reproduce below 
the formulation of LA-SUD: 

max WkXk (5a) 
s. t. — > 7fc \/kE IC, (5b) 

/ J Pm'^mk-^m ~r ^ 
m^k 

Xk E {0, 1} VA; E IC. (5c) 

The objective function (l5al) aims to maximize the total weight of the LA set. The constraints 
(l5bl) formulate the SINR criteria. If = 1, indicating that link k is active, the /cth inequality 
constrains the SINR of link k to be at least 7^. For the case that link k is not active, Xk = 0, 
the A;th inequality in (l5bl) is always satisfied, i.e., it has null effect, if parameter Mk is set 
to a sufficiently large value. By construction, an obvious choice is Mk = Yl PmGmklk + 

m^k 

Vlk — PkGkk- Note that the size of the formulation ([5]), both in the numbers of variables and 
constraints, is of 0{K). 

Now, consider the same system but with receivers having MUD capability that enable 
cancellation of strongly interfering links. We distinguish between IC in a single stage (PIC) 
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and in multiple stages (SIC). In the former, to cancel the transmission of an interfering link, 
all other signals of active links, including the signal of interest, are considered to be additive 
noise, independent of other cancellation decisions at the same receiver. A formal definition 
of the output is given below. 

Problem LA-PIC: Optimal link activation with parallel interference cancellation. 
Output: An activation set ^ C /C and the set ^ A\{k} of cancelled transmissions for 
each k e A, maximizing ^ Wk and satisfying the conditions: 

> 7m Vm e Ck, VA; e A, (6a) 



E PnGnk + V 
nGA\{m} 



>7fc VfceA (6b) 



PmGmk + ri 
m£A\{k,Ck} 



The set of conditions (|6al) ensures that the specified cancellations can take place. For the 
receiver of link k to cancel the transmission of link m, the receiver of k acts as if it was 
the receiver of m. Hence, the "interference-to-other-signals-and-noise" ratio incorporates in 
the numerator the received power PmGmk of the interfering link m and in the denominator 
the received power PkGkk of own link k. This ratio must satisfy the SINR threshold of the 
signal m to be decoded. The set of conditions (l6bl) formulates the SINR requirements for 
the signals of interest, taking into account the effect of IC in the SINR ratio. That is, the 
cancelled terms are removed from the sum in the denominator, determining the aggregate 
power of the undecoded interference which is treated as additive noise. 

For SIC, the output must be augmented in order to specify, in addition to the cancellations 
Cfc, by the receiver of A; G ^, the order in which they take place. A formal definition of the 
output is given below. 

Problem LA-SIC: Optimal link activation with successive interference cancellation. 
Output: An activation set ^ C /C and the set Ck ^ A \ {k} of cancelled transmissions along 
with a bijection : Ck ^ {1, ■ ■ ■ ,\Ck\] for each k E A, maximizing Yl ^k and satisfying 

k(^A 

the conditions: 

PmGmk 



Yl PnGnk + V 

neA\{m,qeCk: fefc{g)<6fc(m)} 



> 7m Vm G Cfc, V/c G A, (7a) 



>7. ykeA. (7b) 



Yl PmGmk + V 
meA\{k,Ck} 
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Similarly to LA-PIC, the set of conditions (fTal) formulates the requirement for SIC and 
the set (TTbl) the requirement for decoding the signals of interest, taking into account the 
effect of IC in the SINR ratio. In the output, the cancellation sequence for each k in the 
activation set is given by the bijection bk', the bijection defines a unique mapping of the link 
indices in the cancellation set Ck to the IC order numbers in the cancellation sequence. That 
is, bk{m),m G Ck defines the stage at which link m is cancelled by the receiver of link 
k. The bijection is used in the IC conditions (fTal) . in order to exclude from the sum in the 
denominator, the interference terms that have been cancelled in stages prior to m. 

III. Complexity 

The baseline problem, LA-SUD, is known to be NP-hard; see, e.g., [fTSll . For a combinato- 
rial optimization problem, introducing new elements to the problem structure may change the 
complexity level, potentially making the problem easier to solve. Hence, without additional 
investigation, the NP-hardness of LA-SUD does not carry over to LA with IC. In this section, 
we provide the theoretical result that problems LA-PIC and LA-SIC remain NP-hard, using 
a unified proof applicable to both cases. 

Theorem 1: Problem LA-PIC is NP-hard. 

Proof: We provide a reduction from LA-SUD to LA-PIC. Considering an arbitrary 
instance of LA-SUD, we construct an instance of LA-PIC as follows. For each link k E JC, 
we go through all other links in /C \ {k} one by one. Let m be the link under consideration. 
The power of link k is set to 

/A \ I VmGmk \ i ^ 1 ,c,^ 

■Pu = max I pfc, I _ ^ -Vj/^kkj, (8) 

where e is a small positive constant. By ([8]), the power of k is either kept as before, or 

grows by an amount such that — — — < 7^. Therefore, link k is not able to decode the 

Pk^kk + V 

signal of m, i.e., the IC condition of LA-PIC cannot be satisfied, even in the most favorable 
scenario that all other links, apart from m and k, are inactive. 

After any power increase of link k, we make sure that this update does not have any effect 
in the application of ([8]) to the other links. This is achieved by scaling down the channel gain 

Gkm as 

G'km - GkmPk/p'k, (9) 

meaning that for any m, the received signal strength from k remains the same as in the 
original instance of LA-SUD. As a result, even though IC is allowed in the instance of LA- 
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PIC, no cancellation will actually take place, since none of the IC conditions holds due to 
the scalings in ([8]) and (|9l). 

By the construction above, for each link k E IC the total interference that is treated as noise 
in the instance of LA-PIC equals that in the instance of LA-SUD. Thus, the denominator of 
the SINR of the signal of interest does not change. On the other hand, the numerator may 
have grown from to p^. To account for this growth, the SINR threshold 7^ is set to 



In effect, the increase of the power on a link, if any, is compensated for by the new SINR 
threshold. Note that, because p^Pk > 1, 7^ prohibits cancellation of the kth signal by any 
receiver other than the kth one. 

From the construction, one can conclude that a LA set is feasible in the instances of LA- 
SUD, if and only if this is the case in the instance of LA-PIC. In addition, the reduction is 
clearly polynomial. Hence the conclusion. ■ 

Corollary 1: Problem LA-SIC is NP-hard. 
Proof: The result follows immediately from the fact that, in the proof of Theorem [H 
the construction does not impose any restriction on the number of links to be cancelled, nor 
to the order in which the cancellations take place. ■ 

IV. Link Activation with Parallel Interference Cancellation 

In this section, we propose a compact ILP formulation for LA-PIC. In addition to the 
Xfc, V/c G /C, variables in (|5]), we introduce a second set of binary variables, Umk, Vm, k E 
/C, m ^ k. Variable ymk is one if the receiver of link k decodes and cancels the interference 
from link m and zero otherwise. The output of LA-PIC is then defined hy A = {k E K, : 
Xk = 1} and Ck = {m E A \ {k} : ymk = 1}, for each k E A. The proposed formulation 
for LA-PIC is 



I'k = IkPk/p'k- 



(10) 



max 




(11a) 



s. t 



Vrnk — 



m 



Vm, k E IC, m ^ k, 



(lib) 



Umk — -^k 



\fm, k E IC, m ^ k, 



(11c) 



PkGkk + Mfc(l - Xk) 



> Ik 



yk E IC, 



(Ud) 



5^ PmCrmki^'^m Vrnk) ~^ V 
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PmGmk + Mmk{\ — ymk) 
PnGnkXn + V 



>7, 



m 



\/m, A; G /C, m ^ k, 



(lie) 



Vmk e {0, 1} 



Vm, A; G /C, m ^ k, 



(llf) 



Xfc G {0,1} 



VA; G /C. 



(llg) 



The objective function (II lal) is the same as (l5al) for LA-SUD. The first two sets of 
inequalities, (II Ibl) and (II Icl) . pose necessary conditions on the relation between the variable 
values. Namely, a cancellation can take place, i.e., ymk = 1, only if both links k and m are 
active, i.e., = = 1- The set of inequalities (llldl) formulates the SINR requirements 
(|6bl) for decoding the signals of interest, in a way similar to (|5bl) for LA-SUD, with the 
difference that here the cancelled interference terms are subtracted from the denominator 
using the term Xm — y-mk- Note that, without (II Ibl) . the formulation will fail, as in (II Idl) it 
would allow to reduce the denominator of the ratio by subtracting non-existing interference 
from non-active links. The next set of constraints (lllel) formulates the condition (l6al) for 
PIC: ymk can be set to one only if the interference from link m, PmGmk, is strong enough 
in relation to all other active signals, including the signal of interest. If the ratio meets the 
SINR threshold 7^ for link m, cancellation can be carried out. Setting ymk to zero is always 
feasible, on the other hand, provided that the parameter Mmk is large enough. A sufficiently 
large value is Mmk — Yl PnGnklm + Vim — PmGmk- The construction of (II lei) reflects the 



fact of performing all cancellations in a single stage, as in cancelling the signal of link m, 
other transmissions being cancelled in parallel are treated as additive noise. Note that the 
model remains in fact valid even if (lllcl) is removed. Doing so would allow the receiver of 
an inactive link to perform IC. However, since an inactive link does not contribute at all to 
the objective function, this is a minor "semantic" mismatch that can be simply alleviated by 
post-processing. 

For practical purposes, each receiver may be allowed to cancel the signal of at most one 
interfering link. The resulting LA problem, denoted LA-SLIC, can be easily formulated by 
adapting the formulation (fTTI) for LA-PIC. The only required change is the addition of the 
set of constraints 




(12) 



m^k 



that restricts each receiver to cancel at most one interfering transmission. 
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The size of the formulation (fTTI) . both in the numbers of variables and constraints, is of 
0{K'^). Thus, the formulation is compact and its size grows by one magnitude in compar- 
ison to ©. In fact, to incorporate cancellation between link pairs, one cannot expect any 
optimization formulation of smaller size. 

When implementing the formulation, two pre-processing steps can be applied to reduce the 
size of the problem and hence speed-up the calculation of the solution. First, the links that are 
infeasible, taking into account only the receiver noise, are identified by checking for every 
receiver whether the received SNR meets the SINR threshold for activation. If the answer 

is "false", i.e., < 7^, then link k is removed from consideration by fixing the Xk 

V 

variable to zero. Second, the link pairs for which cancellation can never take place are found 
by checking, for every receiver and interfering signal, whether the "interference-to-signal- 
of-interest-and-noise" ratio meets the SINR threshold for decoding the interference signal. If 
the answer is "false", i.e., ^J"^™^ < 7^, then link k cannot decode the interference from 

PkCrkk + V 

m and this option is eliminated from the formulation by setting the respective variable ymk 
to zero. 

V. Link Activation with Successive Interference Cancellation Under 

Common SINR Threshold 

Incorporating the optimal IC scheme, SIC, to the LA problem is highly desired, since it 
may activate sets that are infeasible by PIC. However, using ILP to formulate compactly the 
solution space of LA-SIC, is challenging. This is because the formulation has to deal, for 
each link, with a bijection giving the cancellation sequence. We propose an ILP approach 
and present it in two steps. In this section, we consider LA-SIC under the assumption that 
all links have a common SINR threshold for activation, i.e., 'jk = 7, V/c G /C. In the next 
section, we address the general case of individual SINR thresholds. 

For SIC under common SINR threshold, we exploit the problem structure and show that 
the optimal cancellation order can be handled implicitly in the optimization formulation. As 
a result, we show that LA-SIC can in fact be formulated as compactly as LA-PIC, i.e., using 
0{K'^) variables and constraints. The idea is to formulate an optimality condition on the 
ordering of IC. To this end, consider an arbitrary link k and observe that meeting the SINR 
threshold for decoding the signal of interest is equivalent to having in the receiver, after 
IC, a total amount of undecoded interference and noise at most equal to PkGkk/l- We refer 
to this term as the interference margin Uk ■ Similarly, the interference margin that allows 
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cancellation of the interference from link m at the receiver of link k is Umk — PmGmkh- 
Consider any two interfering links m and n, and suppose Umk > Unk- Note that, because the 
SINR threshold is common, the condition is equivalent to PmGmk > PnGnk, i-e-, the receiver 
of link k experiences stronger interference from m. If the condition holds, the cancellation 
of m should be "easier" in some sense. Thus, one may expect that if k can decode both 
m and n, the decoding of m should take place first. In the following, we prove a theorem, 
stating that this is indeed the case at the optimum — there exists an optimal solution having 
the structure in which if a weaker interference signal can be cancelled, then any stronger one 
is cancelled before it. 

Theorem 2: If Umk > Unk and the receiver of link k is able to cancel the signal of n, then 
it is feasible to cancel the signal of m before n and there exists at least one optimum having 
this structure in the cancellation sequence. 

Proof: Let Ink denote the total power of undecoded interference and noise when the 
receiver of k decodes the signal from n. Assume that m has not been cancelled in a previous 
stage. Then, PmGmk is part of Ink- Successful cancellation of n means that Ink < Unk- Since 
Unk < Umk, it holds that Ink < Umk- Considcr now decoding the signal of m immediately 
before n. Thus for this cancellation, the total power of the undecoded interference and noise 
incorporates the interference of n, but not that of m, i.e., Imk = Ink + PnGnk — PmGmk- 
Because Umk > Unk implies PmGmk > PnGnk, it holds that Imk < Ink- Since Ink < Umk, 
the cancellation condition I^k < Umk is satisfied. After cancelling m, IC can still take place 
for n, because the new is decreased by PmGmk- Consequently, both m and n can be 
cancelled. Obviously, doing so will not reduce the number of active links and the theorem 
follows. ■ 

By Theorem [2l for each link k, we can perform a pre-ordering of all other links in 
descending order of their interference margins. SIC at link k can be restricted to this order 
without loss of optimality. At the optimum, the cancellations performed by k, for interfering 
links that are active, will follow the order, until no more additional cancellations can take 
place. In this optimal solution, when considering the cancellation condition of interfering 
link m, interference can only originate from links appearing after m in the sorted sequence. 

We propose an ILP formulation based on Theorem [21 The formulation uses the same 
variables of (fTTl) for LA-PIC, as there is no need to formulate the cancellation order explicitly. 
The sorted sequence is denoted by, for each link G /C, a bijection : JC\{k} t—^ {1, . . . , K— 
1}, where ifc(m) is the position of link m in the sorted sequence. The sorting results in 
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ikim) > ik{n) if Umk < Unk- In case of Umk = Unk, the tie can be broken arbitrarily without 
affecting the optimization result. In addition, let Cmk = K — 1 — ikirn) denote the number of 
links appearing after m in the sorted sequence for k. The proposed formulation for LA-SIC 
under common SINR threshold is 



max ^^WfcXfc (13a) 

keic 

s. t. Umk < Xm Vm, k E }C, m k, (13b) 

Vrak < Xk Vm, k e IC, m k, (13c) 

PkGkk + Mk{l-Xk) ^^^^^ ^^3^^ 



Yl PmGmk{Xm " Umk) + V 
m^k 

PmGmk ~\~ -^rakiX Vrnk) 

PnGnkXn + PkGkk + V 

n^k, i,,{n)>ik(m) 



>7 Vm, /c G /C, m ^ k, (13e) 



^ Vnk < Cmki^ - X-m + Vrnk) ^171, k G JC, Ul ^ k, (13f) 

n^k, ik(n)>ik{m) 

Vrnk e {0, 1} \fm,kelC, m ^ k, (13g) 

Xk e {0, 1} Vfc e /C. (13h) 

The first three constraint sets (I13bl) - (ll3dl) have the same meaning with (lllbh - dlldl) for LA- 
PIC; see Section HYI The constraint set (I13el) formulates the conditions (TTal) for SIC, making 
use of Theorem [21 Consider the condition for cancellation of signal m from receiver k in stage 
ikirn). Then, in the denominator of the ratio, the sum of undecoded interference is limited to 
the transmissions coming after m in the sorted sequence of k, since all other active links with 
higher interference margin than m have already been cancelled. The formulation is however 
not complete without (I13fl) . This set of constraints, in fact, ensures the optimality condition 
set by Theorem [2] and utilized in (I13el) . That is, if both m and n are active, Umk > Unk, and 
n is cancelled by k, then m is cancelled by k as well. Equivalently speaking, if m is active 
but not cancelled by k, then none of the other links after m in the sequence of k may be 
cancelled. Examining (I13fl) . we see that it has no effect as long as Xm equals ymk- If link 
m is active but not cancelled, corresponding to = 1 and y^k = 0, the right-hand side 
of (I13fl) becomes zero, and therefore no cancellation will occur for any n having position 
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after m in the ordered sequence. Also, note that the case Xm = but ymk = 1 cannot occur, 
because of (I13bl) . 

Given a solution to the formulation (fTSl) . the cancellation sequence of each active link 
k, i.e., the bijection bk in the definition of LA-SIC in Section [III is easily obtained by 
retrieving from the predefined bijection 4 the elements with y^m = 1- The compactness of 
the formulation (fT3]) is manifested by the fact that its size, in both the numbers of variables 
and constraints, is of 0{K'^). Thus, provided that there is a common SINR threshold for 
activation, we have formulated LA-SIC as compactly as LA-PIC. 

When implementing the formulation (fT3]) . similar pre-processing steps with (fTTI) for LA- 
PIC can be applied to reduce the size of the problem. First, the infeasible links are removed 
for consideration by fixing Xk to zero when Uk < rj. Second, the infeasible IC options are 
eliminated from the formulation by fixing y^k to zero when Umk < PkGkk + V- 

VI. Link Activation with Successive Interference Cancellation Under 

Individual SINR Thresholds 

In this section, we consider the LA-SIC problem under the most general setup; namely 
when the links have individual SINR thresholds. Differently from the common SINR case, 
treated in Section |Vl a pre-ordering of the sequence of potential IC does not apply. The 
reason is that the interference margin Umk does not depend anymore only on the received 
power PmGmk but also on the link-specific SINR threshold 7^. To see this point, consider a 
scenario where link k attempts to cancel the signal of two interfering links m and n in two 
consecutive stages. Denote by / the sum of the remaining interference, other than m or n, the 
received power of link A;'s own signal, and noise. Assume a mismatch between the relation 
of interference margin and that of received power: PmGmk < PnGnk but Umk > Unk because 
7m < 7n- If cancels m and then n, the cancellation conditions are Umk > PnGnk + I 
and Unk > I- Reversing the cancellation order leads to the conditions Unk > PmGmk + I 
and Umk > I- For our example, we let / = 0.5. Consider two sets of values for the other 
parameters. The values in the first set are: Umk = 3, Unk = 1, PmGmk = 1, PnGnk = 2 and 
in the second set are: Umk = 2, Unk = 1j PmGmk = 0.5 PnGnk = 2. For the first set, both 
interfering links can be cancelled only if cancellation applies to m first, whereas the opposite 
order must be used for the second set. Hence the interference margin (or received power) 
does not provide a pre-ordering for cancellation. 

In the following, we propose an ILP formulation for LA-SIC under individual SINR 
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thresholds, that explicitly accounts for the cancellation order. Our approach is to introduce 
for each pair of links, m, k E IC, a set of binary variables y^^, and represent the cancellation 
stage by the superscript t. Variable y^;. is one if and only if the receiver of link k cancels the 
interference from link m at stage t. The effect is that, for each link k, the solution values of 
yl^i^ order the feasible cancellations; hence, they have a direct correspondence to the output 
bijection bk of LA-SIC, defined in Section [III It is apparent that the index t ranges between 
one and K — 1. In practice, due to computational considerations, we may want to restrict 
the maximum number of cancellation stages. To this end, we define, for each k E IC, the 
integer parameters Tk < K — 1 and the sets Tk = {I, . . . ,Tk}. The proposed formulation of 
the general LA-SIC problem, under individual SINR thresholds and restricted cancellation 
stages, is 



max ^^WfcXfc (14a) 

S. t. 

^ Vmk < a^m, Vm, k eK:, m^k, (14b) 
t=i 

5^ VL < yk G /C, Vt G Tk, (14c) 

PkGkk + Mk{l-Xk) ^^^^^ ^^^^^ 



E Pm.Gmk[Xm -YvLk) +V 
mj^k ^ t=l ' 



PmGmk + MmkiX " 



mk) 



Yl PnGnk [Xn- Y Vtk ) + PkGkk + Tj 
n^m,k ^ t'=l ' 



> 7m Vm, k E }C, k, ETk, (14e) 



J2yLk<J2y'^^' v/c G /c, vt G ra {1}, (m 

mj^k mj^k 

ylk^{0,l} ym,kEJC, m^k,ytETk, (14g) 

Xk E {0, 1} VA; G IC. (14h) 

The conditions (I14bl) - (ll4cl) have similar role with (I13bl) - (ll3cl) . Namely, only when links 
k and m are active, the receiver of k can consider to cancel the transmission of m. In 
addition, the summation over t in (I14bl) ensures that link m is cancelled in at most one 
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Stage. Furthermore, the summation over m in (I14cl) enforces each receiver to perform at 
most one cancellation per stage. This removes equivalent solutions, without compromizing 
optimality, to enhance computational efficiency. The SINR requirements for decoding the 
signals of interest are set in (I14dl) . similarly to (I13dl) . In the denominator, all cancelled links, 
regardless of the stage the cancellation is performed, are removed from the sum of undecoded 
interference. The next set of constraints (I14el) formulates the requirement for the cancellation 
of link m by link k at stage t. The active interfering transmissions that have been cancelled 
before stage t are excluded from the sum of undecoded interference in the denominator of 
the ratio. The constraints (I14el) are formulated with the convention that the sum within the 
parenthesis in the denominator of the ratio is zero for t = 1. Note that, even though for 
each receiver k and interfering link m, constraints are formulated, due to (I14bl) . all but 
at most one will be trivially satisfied by the respective variabls being equal to zero. 
The constraints (I14fl) are not mandatory for the correctness of the formulation, but their role 
is to enhance the computational efficiency. These constraints ensure that the cancellations 
are performed as "early" as possible, i.e., there are no "idle" stages at which no cancellation 
takes place and which are followed by later stages where cancellation takes place. Otherwise, 
if mi and m2 are cancelled by k, and the cancellation of the former takes place first, the 
cancellations can be performed at any two stages ti and t2, as long as ti < ^2- Clearly, 
such solutions are all equivalent to each other, and the presence of them slows down the 
computational process. 



Since (fT4l) formulates the most general LA-SIC problem, it also applies to the common- 
SINR case of Section |Vl Its computational efficiency though is significantly lower than the 
respective of formulation (fT3]) . The reason is that its size is one magnitude larger than (fT3]) . 
i.e., the numbers of variables and constraints grow from 0{K'^) to 0{K^). However, we note 
that the formulation (fT4l) remains compact. In order to deal with the scalability issue, one may 
resort to restrict the maximum number of cancellations, T^., to a constant being considerably 
lower than K—1. Typically, doing so has little impact on the solution quality, because most of 
the performance gain from IC is due to the first few cancellations. Also, when implementing 
the formulation (fT4l) . similar pre-processing steps with ([13]) can be applied, see Section |Vl 
to reduce the size of the problem. 
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(a) I dataset; sparse topology 



(b) N dataset; sparse topology 
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(c) I dataset; dense topology 



(d) N dataset; dense topology 



Fig. 1. Instances of a 20-link network for different datasets and densities; transmitters marked with circles, receivers with 
triangles. 



VII. Numerical results 

This section presents a quantitative study of the effect that IC has on the optimal LA 
problem in wireless networking. The ILP formulations, proposed in Sections HVl - IVIl are 
utilized to conduct extensive simulation experiments on randomly generated network instances 
with various topologies, densities, cardinalities, and SINR thresholds. Nodes are uniformly 
scattered in square areas of 1000m x 1000m and 500m x 500m, in order to create sparse 
and dense topologies, respectively. Two types of datasets are generated. The first one takes 
an information-theoretic viewpoint and is henceforth denoted dataset I. To this end, the 
transmitter-receiver matchings are arbitrarily chosen |[8l, EOl . with the sole criterion of 
feasible single-link activation. Thus, the links have arbitrary length within the test area, 
provided that their SNR is larger than the SINR threshold required for activation. The second 
dataset provides a rather networking-oriented approach [fTSll . OTll and is henceforth denoted 
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I dataset; sparse topology 
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(a) I dataset; sparse topology 

I dataset; dense topology 
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N dataset; sparse topology 
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(b) N dataset; sparse topology 

N dataset; dense topology 
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(c) I dataset; dense topology 
Fig. 2. Distribution of SNR for different datasets and densities. 



SNR (dB) 

(d) N dataset; dense topology 



dataset N. In this dataset, the length of the links is constrained to be from 3m up to 200m, with 
the rationale to produce instances resembling a multihop network. The networks considered 
have cardinality K ranging from 5 up to 30 links. Fig. [T] illustrates instances of a 20-link 
network; Figs. [Tal and [Tel correspond to the sparse and dense topology, respectively, of dataset 
I, whereas Fi g s . [Tbl and [Tdl correspond to the sparse and dense topology, respectively, of dataset 
N. 

The input parameters are chosen to be common for all links; specifically, the transmit 
power pk, VA; G /C, is set to 30dBm, the noise power i] to — lOOdBm, and the channel 
gains G„ik follow the geometric, distance-based, path loss model with an exponent of 4. The 
major difference between the datasets is the distribution of the link lengths, which effectively 
determines the SNR distribution of the links. The input parameters yield minimum SNR 
approximately equal to 4dB, 16dB, and 32dB for dataset sparse I, dense I, and N, respectively. 
The histograms in Fig. [2] illustrate the SNR distribution of each dataset; as in Fig. [H left and 
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(c) LA-PIC (d) LA-SIC 

Fig. 3. Exemplary activation sets for different IC schemes; I dataset; dense topology; 10 links; SINR — 9dB. 



right sub-figures are for dataset I and N, respectively, whereas upper and lower sub-figures are 
for sparse and dense topologies, respectively. For dataset I, the links in the sparse topology 
have on average lower SNR than in the dense topology; the mass of the SNR distribution is 
roughly for 10-40dB in the sparse and for 20-50dB in the dense topology. This is because 
in the sparse topology the test area is enlarged, allowing generation of longer links which 
have lower SNR values. On the contrary, the SNR distribution of dataset N is invariant to 
the network density; this is by construction, since the distribution of the link lengths is not 
affected by the size of the test area. 

For each dataset and network cardinality, 30 instances are generated and the performance 
of LA with IC is assessed by two simulation studies. In the first study, all links are assumed 
to require for activation a common SINR threshold 7^ = 7, V/c G /C, taking values from 
— 9dB up to 6dB, and have equal activation weights, e.g., Wk = 1, VA; G IC. The goal is 
to evaluate the performance gain due to single-link, parallel, and successive IC schemes on 
the LA problem over the baseline approach without IC. For this purpose, we implemented 
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SINR -6 dB; I dataset; sparse topology; 30 instances 



30 




_o_ 


LA- 


-SUD 






LA- 


-SLIC 


« 25 
c 


- * - 


LA- 
LA- 


-PIC 
-SIC 


■a 

0) 








1 20 








(0 








15 








03 
SI 

E 
















S 5 


n — 




--S 











5 




10 



SINR -6 dB; N dataset; sparse topology: 30 instances 
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SINR -6 dB; I dataset; dense topology; 30 instances 
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(b) N dataset; sparse topology 
SINR -6 dB: N dataset; dense topology: 30 instances 
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Fig. 4. Average number of activated links versus network size for SINR threshold — 6dB. 



the formulations ©, (fTT])-(fT2l). and for LA-SUD, LA-SLIC, LA-PIC, and LA- 
SIC, respectively. Fig. |3] illustrates exemplary activation sets for an instance of a 10-link 
network, drawn from dataset dense I, when the SINR threshold is — 9dB. It is evidenced that 
performance increases with problem sophistication: Figs. (I3al) . (I3bl) . (I3cl) . and (l3dl) . show that 
LA-SUD, LA-SLIC, LA-PIC, and LA-SIC activate 3, 5, 6, and 10 links, respectively. 

The optimal solutions are found by an off-the-shelf solver, implementing standard tech- 
niques such as branch-and-bound and cutting planes [[5l. The simulations were performed 
on a server with a quad-core AMD Opteron processor at 2.6 GHz and 7 GB of RAM. 
The ILP formulations were implemented in AMPL 10.1 using the Gurobi Optimizer ver. 
3.0. Regarding the computational complexity of the proposed ILP formulations for IC, an 
empirical measure is the running time of the solution process. We have observed that it is 
not an obstacle for practical instance sizes. 
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SINR 3 dB: I dataset; sparse topology: 30 instances SINR 3 dB; N dataset; sparse topology; 30 Instances 
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(a) I dataset; sparse topology (b) N dataset; sparse topology 

Fig. 5. Average number of activated links versus network size for SINR threshold 3dB. 



In the following a selection of the simulation results is presented. Fig. |4] shows the average, 
over 30 instances, number of activated links versus the total number of links in the network, 
achieved by all versions of the LA problem when the SINR threshold is — 6dB. The results in 
the four sub-figures correspond to the datasets exemplified in Fig. \T\ The major observation 
is that all LA schemes with IC clearly outperform LA-SUD and in particular LA-SIC yields 
impressive performance. Comparing Figs. |4a] and |4cl it is concurred that the results for dataset 
I are density invariant. As the number of links in the network increases, the performance of 
LA-SUD improves, due to the diversity, almost linearly but with very small slope. LA-SIC 
though improves significantly, activating two to three times more links than the baseline. 
When the network has up to about 15 links, nearly all of them are activated with LA-SIC. 
On the other hand, LA-PIC has a consistent absolute gain over LA-SUD, activating one to 
two links more. Furthermore, LA-SLIC has almost as good performance with LA-PIC, i.e., 
it captures most of the gain due to single-stage IC. Fig. |4d] shows that the LA schemes have 
similar performance in dataset dense N as in dataset I. Fig. |4b] shows that LA is easier for 
dataset sparse N, even without IC. The curves of all schemes linearly increase with network 
cardinality, but with IC the slopes are higher, so that the absolute gains, differences from 
the baseline, broaden. Maximum gains are for 30 links, where LA-SUD, LA-SLIC, LA-PIC, 
and LA-SIC activate about 18, 22, 23, and 30 links, respectively. For the tested network 
cardinalities, LA-SIC achieves the ultimate performance activating all links. 

As seen in Fig. H the performance gains due to IC are very significant when the SINR 
threshold for activation is low. However, for high SINR thresholds, the gains are less promi- 
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Fig. 6. Average number of activated links versus SINR threshold for network of 30 links. 



nent. For example, Fig. [5] shows the performance of the LA schemes when 7 is set to 3dB. 
Figs. [5a] and [5b] are for the sparse datasets I and N, respectively; for the dense topologies the 
results are similar to Fig. [Sa] It is evidenced that IC schemes activate one to two links more 
than the baseline and that most of this gain can be achieved with single-stage IC. 

The fact that the IC gains diminish with increasing the SINR threshold is clearly illustrated 
in Fig. [6l which compares, for networks of 30 links, the average performance of all LA 
schemes for various SINR thresholds. The relative gain of SIC is more prominent in the case 
of dataset I, which is more challenging for the baseline problem. For dataset I, when the 
SINR threshold is low, around — 9dB, SIC activates nearly all links, whereas SUD activates 
less than a third of them. For sparse and dense dataset N, SIC activates effectively all links 
when the SINR threshold is lower than — 6dB and — 9dB, respectively, whereas SUD activates 
less than two thirds and less than half of them, respectively. For mid-range SINR thresholds, 

March 17, 2012 DRAFT 



24 



Submitted to IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY on March 17, 2012 



20 
„ 18 

N 
I 

« 16 

£^ 

Z 10 

o 

I 8 
c 

03 E 
Ul t> 

2 



I dataset; sparse topology; 30 instances 



-no IC 
-1 IC 
-2 ICS 
-3 ICS 
-4 ICS 
-5 ICS 




15 20 25 

number of available links 



N dataset; sparse topology; 30 instances 



-no IC 
-1 IC 
-2 ICS 
-3 ICS 
-4 ICS 
-5 ICS 




(a) I dataset; sparse topology 



10 15 20 25 

number of available links 

(b) N dataset; sparse topology 



20 
„ 18 

N 
X 

3 16 
t 10 



I dataset; dense topology; 30 instances 



-no IC 
-1 IC 
-2 ICS 
-3 ICS 

-4 ICS 
-5 ICS 




10 15 20 25 

number of available links 

(c) I dataset; dense topology 



N dataset: dense topology; 30 instances 




10 15 20 25 

number of available links 

(d) N dataset; dense topology 



Fig. 7. Average network throughput versus network cardinality for various IC schemes and SINR thresholds randomly 
chosen from {—6, —3, 3}dB. 



up to about 3dB, SIC has an exponentially decreasing performance, but nevertheless still 
significantly outperforms SUD. On the other hand, for SINR thresholds up to about OdB, PIC 
yields a relatively constant performance improvement of roughly two to five links, depending 
on the dataset. PIC is effectively equivalent to its simpler counterpart SLIC, for SINR higher 
than — 6dB. The performance of all IC schemes converges for SINR thresholds higher than 
3dB. The interpretation is that if IC is possible, it is more likely that it will be restricted to 
a single link. For very high SINR thresholds, it becomes rarely possible to perform IC. 

In the second simulation setup, the performance of the general LA-SIC problem, under 
individual SINR thresholds, is evaluated. The SINR threshold 7^ for each link is taking, with 
equal probability, one of the values in the set {—6, — 3,3}dB and the activation weights Wk 
are set equal to the data rates, in bits per second per Hertz, corresponding to the respective 
SINR thresholds. The formulation (fl?)) is implemented varying the maximum number of 
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cancellation stages Tk = T , \/k ^ /C, from 0, corresponding to the baseline case without 
IC, up to 5. Fig. |7] shows the average, over 30 instances, throughput of all activated links 
versus the network cardinality, for all the datasets. For dataset I, the network throughput is 
almost doubled with IC; roughly half of this increase is achieved by the first cancellation 
stage and most of the rest by the next two to three stages. For dataset N, it is seen that the 
first cancellation stage yields a significant gain of about 2 b/s/Hz and that it only pays off 
to have more than two cancellation stages for large and sparse networks. 

VIII. Conclusions 

In this paper, we have addressed the problem of optimal concurrent link activation in 
wireless systems with interference cancellation. We have proved the NP-hardness of this 
problem and developed integer linear programming formulations that can be used to ap- 
proach the exact optimum for parallel and successive interference cancellation. Using these 
formulations, we have performed numerical experiments to quantitatively evaluate the gain 
due to interference cancellation. The simulation results indicate that for low to medium 
SINR thresholds, interference cancellation delivers a significant performance improvement. 
In particular, the optimal SIC scheme can double or even triple the number of activated 
links. Moreover, node density may also affects performance gains, as evidenced in one of 
the datasets. Given these gains and the proven computational complexity of the problem, 
the development of approximation algorithms or distributed solutions incorporating IC are of 
high relevance. 

Concluding, the novel problem setting of optimal link activation with interference can- 
cellation we have introduced here provides new insights for system and protocol design in 
the wireless networking domain, as in this new context, strong interference is helpful rather 
than harmful. Thus, the topic calls for additional research on resource allocation schemes in 
scheduling and routing that can take the advantage of the interference cancellation capability. 
Indeed, the LA setup studied herein assumes fixed transmit power for active links. This can 
lead to increased interference levels, since the SINRs can be oversatisfied. Incorporating 
power control to the LA problem with IC will bring another design dimension which can 
yield additional gains. Furthermore, it may enable IC for high SINR thresholds IC. 
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